On the theory of BoseEinstein condensation of 
quasiparticles: on the possibility of condensation of 
ferromagnons at high temperatures 

A.I. Bugrij 1 and V.M. Loktev 2 

N.N. Bogolyubov Institute for Theoretical Physics 
of the National Academy of Sciences of Ukraine 
ul. Metrologicheskaya lJ^-b, Kiev 03143, Ukraine 
E-mail: 1 abugrij@bitp. kiev.ua, 2 vloktev@bitp. kiev.ua 

Fiz. Nizk. Temp. 33, 5168 (January 2007) 
(Submitted July 31, 2006) 

Abstract 

The Bose condensation of magnons in physical systems of finite size is consid- 
ered for the case of ferromagnetic thin films. It is shown that in accordance with 
present-day experimental capabilities, which permit one to achieve densities of long- 
wavelength spin excitations of ~ 10 18 - 10 19 ~ 3 , in such films, the formation of a 
coherent condensate of such quasiparticles begins at temperatures T ~ 10 2 K (in- 
cluding room temperature). It is found that Bose condensation is accompanied 
by a scaling phenomenon, according to which the main thermodynamic variable is 
not the number of particles N but the ratio N/T . This indicates that the Bose 
condensation of magnons can be observed at relatively low magnon densities (and, 
accordingly, low pumping). The roles played by the shape of the spectrum of spin 
excitations and by the film thickness for observation of the phase transition to the 
state with the Bose condensate are analyzed, and the partial contributions of dierent 
groups of quasiparticles to the total spectral distribution of magnons over energies 
are discussed. 

PACS: 05.30.Jp, 75.30.Ds, 75.70-i 
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1 Introduction 



The BoseEinstein condensation (BEC) of atoms and molecules pfl [2] has become one of the 
most remarkable phenomena that reveal and confirm the quantum nature of a number 
of macroscopic processes. The formation of a Bose condensate, i.e., the accumulation 
of identical particles with integer spin in one of the quantum states, may be inherent 
both to true particles (atoms, molecules) and to quasiparticle excitations of multiparticle 
systems. In this sense, quasiparticles — excitons and biexcitons and also magnons are 
of particular interest, since, existing only as excited states, they are actually absent (if 
one ignores the thermal background) in systems at normal temperature and pressure. 
Ordinarily, except under specially chosen conditions, the equilibrium density of thermal 
quasiparticles decreases with decreasing temperature. Therefore, the study of BEC of 
quasiparticles, or their "nonthermal" accumulation in one of the states, requires, first 
and foremost, the presence of a macroscopic (essentially nonequilibrium) total number of 
quasiparticles. This can be achieved only by using methods of creating and maintaining a 
large number of quasiparticles in condensed systems, at least for a time sufficient for, first, 
their relatively rapid thermalization and, second, subsequent BEC. Under such conditions 
the latter will occur as a (quasi-) equilibrium phenomenon at a conserved (on average) 
number of quasiparticles, which is ensured by some external source of intense quasiparticle 
creation. 

BEC is a phenomenon associated, as a rule, with very low temperatures: the critical 
temperature T B ec & t which a Bose condensate is created depends on the gas density n 
and the mass m of its constituent particles according to the well-known formula [3] 

n 2/3 2-rrh 2 

Tbec = ^^3/2) ' (L1) 

where ((x) is the Riemann ( -function. For example, in experiments on the observation 
of BEC jH E] because of the large mass of the atoms of alkali elements and low density 
(approximately 10 3 particles in a volume of ~ 1CT 6 3 ) the system must be "cooled" 
down to . . . 10~ 8 K. Or, for example, in the case of the "most ideal" of the real 
gases, helium, the temperature of the transition to the Bose-condensed state has a value 
Tbec ~ 2 x 10~ 2 at a density ~ 10 19-3 , which corresponds to a gaseous state. It is 
therefore obvious that observation of BEC at high (up to room) temperatures is possible 
only in systems consisting of light (weakly interacting) bosons. As has been mentioned 
repeatedly, quasiparticle excitations, the effective mass of which is rather small — in 
particular, comparable to the electron mass m e - - are well suited to the role of such 
bosons. 

In the discussion of the BEC of quasiparticles, attention has apparently been focused 
primarily on excitons and biexcitons. Starting with the pioneering work of Moskalenko 
[0] and Keldysh [7], many papers have been devoted to the study of the phase diagram 
of semiconductors at high densities of nonequilibrium electrons and holes and also to 
attempts to observe different types of their joint condensed states. Despite the short 
lifetime of excitonic states and the intense processes of electronhole annihilation, there 
have been several pieces of convincing experimental evidence in favor of the observation 
of exciton (and biexciton) BEC [8]. 

However, there are other promising systems for studying the BEC of quasiparticles. 
In particular, among the interesting objects suitable for this (see Ref. [9]) are magnets, 
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including ferromagnetic insulators, where high densities of magnons (up to 10 18 - 10 19 ~ 3 ) 
can be created by pulsed microwave pumping. We recall that magnons, which are el- 
ementary excitations above the magnetization field, to good approximation obey Bose 
statistics^ Their spectrum is formed on account of the presence of several interactions, 
chief among which are the exchange and magnetic dipoledipole interactions [10]. The first 
of these determines an isotropic spectrum of spins waves with wavelengths k~ x , less than 
the characteristic size L of the system, so that kL ^> 1 . The second, on the contrary, 
operates at wavelengths for which kL < 1 and, furthermore, leads to dependence of the 
magnon spectrum on the direction of k . Analysis of the conditions for BEC of specifically 
these elementary excitations in comparatively thin slabs (microfilms) is the subject of the 
present study, which is directed toward the task of examining the conditions for obser- 
vation and the features of BEC in ferromagnetic insulating films with different shapes 
of the spin-wave spectrum. We consider magnons to be very promising for research on 
the BEC of quasiparticles for at least several reasons. First, they have a relatively long 
lifetime and, as we have said, can have densities reaching ~ 10 19 ~ 3 ; second, there are 
technologies in place for growing very thin (less than < 10 fim thick) films with rather 
perfect structure and for making precision measurements of the magnon spectrum in them 
by optical methods; third, the temperature and magnetic field can be easily varied over 
wide limits, permitting detailed comparison of theory and experiment. 

We note that in recent years a rather large number of papers have appeared (e.g., 
Refs. [H]-[l6j) in which the idea of BEC of magnons has been used for describing phase 
transitions in (predominantly) antiferromagnets from their nonmagnetic (singlet) state to 
a magnetically ordered state under the influence of an external magnetic field. The point 
is that the induced appearance of magnetization in finite fields can actually be formally 
described and interpreted in the language of condensation of magnetic excitations. How- 
ever, BEC as such does not occur, since there one is talking about only a rearrangement 
of the ground state of the system and, consequently, of virtual and not real magnons (see, 
e.g., Ref. [IT]). In contrast, the focus of our attention here is on real — excited — states 
that appear as a consequence of pumping electromagnetic energy into the system, this 
energy going to the creation of quasiparticles above the ground state, as was masterfully 
done in a recent experimental study with the use of microfilms of yttrium iron garnet 

(yig)hbi. 

2 Model and general relations 

Suppose that we have a ferromagnetic crystal in the form a parallelepiped (see Fig. 1) of 
volume V = L x L y L z . The number of sites Nj along the corresponding axes (j takes the 
values x, y, z) is determined by the lattice parameters aj : Nj = Lj/ aj + 1 . As we know, 
the magnetic ground state of such a crystal presupposes an identical direction of all the 
spins (e.g., along the quantization axis z), while the lowest excited state corresponds to 
one "flipped" spin [191 [20]. From a linear combination of such states one can construct an 

1 Interestingly, it was stated in Ref. [H] that all known cases of BEC involve the concept of composite 
bosons formed by an even number of fermions. This pertains to both particles helium and alkalimetal 
atomsand to quasiparticlesCooper pairs, excitons, and biexcitons. This is undoubtedly true if one is 
talking about particles or about WannierMott excitons and biexcitons. It is harder to agree with such a 
statement for small-radius excitons, or Frenkel excitons. And it clearly does not apply to such quasipar- 
ticles as magnons, especially in insulators describable by the Heisenberg model. 
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Figure 1: Shape of the ferromagnetic crystal. 



eigenstate in the form a spin wave with amplitude 6 q (r) . The form of these amplitudes 
depends on the boundary conditions. Usually for simplification one uses cyclic boundary 
conditions, when the amplitudes are identical at opposite faces of the crystal. Then the 
amplitudes have the form of a plane wave: 

6 q (r) = exp(zqr) J]^, (2.1) 

j 

where in the given case the vector r enumerates lattice sites — r 3 - = 1, 2, . . . , Nj , and 
the dimensionless quasimomentum q is defined in the Brillouin zone (— 7r < qj < tt) and 
takes on a discrete spectrum of values with a step Aq^ er = 2n /Nj . 

If, on the contrary, the spins on the faces are "free", the solution of the boundary 
condition problem with the corresponding (free) boundary conditions leads to amplitudes 
that differ from (HTTP : 

6 q (r) = n sin (^')[2/(iV J + l)] 1/2 . (2.2) 

j 

The Brillouin zone here is defined somewhat differently: < qj < tt , and the dis- 
creteness step is Aqj ree = n/(Nj + 1) . Let us point out the main differences between 
the quasimomentum spectra in (12. ip and (12. 2p . First, for periodic boundary conditions 
there exists a zero mode q = , while for the free boundary conditions there is not: 
the minimum value that can be taken on by each of the quasimomentum components 
is qJ Lm = Aq^ Tee . Second, the quasimomentum spectrum for free boundary condi- 
tions is twice as dense as for periodic boundary conditions, since for Nj ^> 1 the ratio 
Agf 7Agf ee = 2N j /{N j + 1) « 2 . 

It is considered to be almost obvious that the concrete form of the boundary conditions 
does not influence the behavior of physical quantities. Actually this is true if one is talking 
about very "large" systems. However, as will be shown below, for "small" systems the 
difference mentioned above — the absence of a zero mode in the spectrum for the case 
of free boundary conditions — leads to quite perceptible contributions (Weyl corrections) 
to certain observable characteristics which are absent in the case of idealized periodic 
boundary conditions. Since we intend to investigate the question of BEC in microfilms 
~ lfim thick, the use of free boundary conditions is justified as more realistic. 
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In spite of the differences indicated, the expression for the magnon energy due to the 
isotropic exchange interaction is the same for both kinds of boundary conditions and can 
be written in dimensionless form as 

e(q) = 2$>n 2 |. (2.3) 

j 

Far from the Brillouin zone boundary or in long-wavelength region (qj <C 7r) , it follows 
from (El) that 

e(q) = y • (2-4) 

The dimensional quantities are easily restored through the dimensions of the crystal 
Lj , the lattice parameters cij , and the effective mass m m , which for magnons is inversely 
proportional to the exchange integral. Then the components of the quasimomentum for 
the case of free boundary conditions 



V . = %^L^ hir^-, kj = 1, 2, . . . , Nj, 



while the dispersion relation 



2m m "i--*> ^ 2m m L] 

3 J 

corresponds to the simplest one for magnons in a ferromagnet. Also taking into account 
the fact that the long-wavelength magnons interact weakly, ~ (P1P2) 2 > with each other 
and (with practically the same amplitude) with phonons|20j. we arrive at the problem of 
BEC of an idea|§ Bose gas, for which the partition function Z m in the grand canonical 
ensemble at temperature T and chemical potential fi m has the form 

In Z m = - H 1 ~ ex Ph(^k - fi m )/T) } , (2.6) 

k 

where 5k — the energy of the state with quantum numbers k = (k x , k y , k z ) , is defined 
in (12.51) . In (12.61) and below we employ the system of units ks = h = 1 , restoring the 
dependence on the fundamental constants as necessary; in addition, the index m will be 
dropped with the understanding that the particles under discussion are magnons. Then 
for the mean number of particles in quantum state k (occupation number) we have 

n k = {ex V [(s k -fi)/T]-l}- 1 , (2.7) 

and for the total (mean) number of particles in the system 

N = T^ = J2n k . (2.8) 

It follows from the domain of definition of the thermodynamic quantities (I2.6l) - (l2.8p 
that the range of variation of the chemical potential for Bose systems is bounded by 



2 The ideality is determined by the density of Bose excitations, which even for a quasiparticle number 
10 20-3 turns out to be extremely small: ~ 10 -3 per atom. 
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the minimum value of the energy in the dispersion relation under study. In the case of 
dispersion relation (12. 5p the energy reaches a minimum value Eq in the quantum state 
with k= (1,1,1): 

71-2 r-2 

3 

We note, by the way, that the widespread assertion that the chemical potential is equal 
to zero as a consequence of the free creation and annihilation of (quasi)particles is not 
completely correct. This error apparently stems from the fact that the free energy F(N, T) 
in the canonical ensemble reaches a minimum as a function of the number of particles (at 
fixed volume and temperature) at \x = . This is essentially a trivial consequence of one of 
the definitions of the chemical potential as a thermodynamic function, [i = OF/ dN . But 
why is the condition of minimum free energy equivalent to the condition of free creation 
and annihilation of particles? Logically it would seem to flow from the requirement of 
maximum entropy as a function of N . Meanwhile, it is easy to check by elementary 
calculations (in the case of an ideal gas, at least) that the entropy is a monotonically 
increasing function of the number of particles and reaches its maximum value at N —>■ oo . 
Analogously, in the grand canonical ensemble the entropy is a monotonically increasing 
function of the chemical potential and reaches its maximum at the maximum admissible 
value of the latter, i.e., when fi — » 6q . Here both the mean number of particles iV and the 
entropy itself go to infinite in a finite volume. The chemical potential /x as an independent 
thermodynamic variable, although a convenient parameter for theorists, is absolutely a 
formal quantity from the standpoint of experimentalists, since there is no prescription 
for its direct measurement. In experiment the chemical potential can only be inferred 
indirectly, e.g., from measurement of the mean number of particles or other observables. 
Then the value of [i can be recovered on the basis of some theoretical prescriptions: as a 
rule, from the formulas of statistical mechanics of an ideal gas. In the case of an ideal gas, 
however, one can eliminate the chemical potential completely from the thermodynamic 
formulas by replacing this independent thermodynamic variable by some other quantity 
having a more direct physical content. If one is talking about BEC, a completely suitable 
candidate for this role [21] is the number of particles at the lowest level, no ■ In fact, it 
follows directly from Eq. (12.71) that 

e (eo-K)/T = 1 + _ (2>g) 

Then, substituting (JOJ into GZM- GZty . we obtain 
InZ = ln(n + l) -J^'ln 

k 

n k = [(l + n -i)e^- e °)/ r -l]- 1 , (2.11) 
N = n + ^'n k , (2.12) 

k 

where the prime on the summation sign means that the term corresponding to the lowest 
magnon state with quantum numbers k = (1, 1, 1) is excluded. 

Such a parametrization is convenient for two reasons. First, it gives a formal definition 
of the condensate in a finite system as simply a set of particles at the lowest-energy 



1 _ e -Ok-eo)/T. 



n 



n + 1 



(2.10) 
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quantum level: the number of these particles is A^bec — n o ■ Furthermore, it permits 
a correct transition to the thermodynamic limit in expressions obtained for systems of 
finite size. This question will be discussed in a little more detail when we specialize to 
microfilms. 

Of course, in the sums over quantum states in (I2.6p . (12.81) and (12.101) . (I2.12p the depen- 
dence on the size of the system enters implicitly through the spectrum of the Schrodinger 
operator, which depends on the boundary conditions. To isolate this dependence we resort 
to the standard technique of changing from summation over quantum states to integration 
over the phase volume: 

£)/(Pk)-> / d$f(p), (2.13) 
k ^ 

where the element of phase volume is usually taken as 

It should not be forgotten, however, that expression (I2.14p is only the first term of an 
asymptotic expansion for p — > oo . According to the famous Weyl problem [22] of the 
number of eigenvalues of an operator which do not exceed a specified value, the coefficients 
of the corresponding asymptotic series are expressed in terms of geometric invariants. In 
the particular case of the Schrodinger operator, taking the next term into account leads 
to the following expression for the phase volume: 

rf$ = d& v - d$ Sv > 
d$ Sv = P^Pdp, (2-15) 

where Sy is the surface area of the sample of volume V . It was mentioned above that the 
contribution d&s v is due precisely to the absence of the zero mode in bounded systems. 
It is not difficult to see that integration over the phase volume elements (I2.14p and (I2.15P 
generates the quasiclassical expansion for the corresponding thermodynamic variables. 
Transformation (I2.13P also presupposes that the function f(p) in the integrand is rather 
smooth. Otherwise some terms (say, the first) can differ sharply in value from the others. 
Then one needs to separate them off and apply transformation (12.131) to the remaining 
sum. A simple example illustrating the difference between the sum and integral in such 
a situation is given in Appendix A. The partition function in the BEC regime, i.e., for 
ji — > Eo , is the case that was indicated in [3], for example. It is seen that the change of 
independent thermodynamic variable (I2.9P from the chemical potential to no solves in 
one stroke the problem of separating out the singular contributions [cf. (12. 6p . (12.81) and 
(I2.10p . (12.121) ]. The corrections to the thermodynamic variables due to d&s v we wm can 
surface contributions, and the contributions due to d§ v , volume contributions. We note 
that their ratio is not universal, in the sense that it can be larger or smaller for different 
thermodynamic functions under equal conditions. In other words, as the volume of the 
system increases, the thermodynamic limit sets in sooner for some physical quantities 
and later for others. The study of rather thin films is clearly a step in the direction of 
the physics of mesoscopic systems, or "nanophysics" , as it has come to be called. As a 
quantitative estimate of the boundary between "macroscopic" and "mesoscopic" one can 
use a comparison of the volume and surface contributions to a given physical quantity. 
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When they are comparable, the formulas of ordinary macroscopic thermodynamics cease 
to work. In particular, the energy, free energy, entropy, and even the mean number of 
particles lose the property of extensivity and, for example, such an interesting quantity 
as the pressure loses the property of isotropy. 

Turning to concrete calculations, we note that averaging over the phase volume (I2.15P 
leads to integrals of the form 

oo 

= r(a)P (0 , (2.16) 

where T(a) is the gamma function, and P a {z) is a polylogarithm, a special function with 
rather simple properties. In the region Re x > it has the series expansion 

oo 

P a (e- x ) = Y. rae ~ lX > ( 2 - 17 ) 

1=1 

which implies that 

4-Pa(e-') = -Pa(*-') ■ (2-18) 
ax 

At x = there is a branch point P a (e~ x ) ~ x a ~ l for noninteger a, and P a (e~ x ) ~ 
x a ~ 1 lnx for integer a; this branch point is explicitly separated out. For example, the 
function P 5 / 2 (e~ x ) has the representation 

JWO = + f; c -^^(-xy , (2.i9) 

1=0 

where the radius of convergence of the series on the right-hand side of where the radius 
of convergence of the series on the right-hand side of (I2.19P is \x\ < 2tt ; as we have said, 
((I) is the Riemann ( function. 

Let us calculate, for example, the mean number of particles (12.131) 

N = n + N ex , 
N cx = ^'n k = N v -N Sv , 

k 

Nv = ^Pw(-^tt), (2-20) 



Sy ( n \ S- 



where 



is the so-called thermal de Broglie wavelength. It follows from (12.201) that at a tem- 
perature of 1 K, quasiparticle mass m e , and number of particles in the condensate 
no ~ 10 16 , the contributions Ny and Ns v become equal when the film thickness de- 
creases to L x ~ 1 fim , which, in accordance with what we have said above, is a direct 
indication of mesoscopicity of the system. 
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The problem of BEC of an ideal gas is one of the few that have an exact solution for 
a large (but finite) number of particles. The solution permits one to trace the formation 
of nonanalyticity of some physical quantity or other as a function of temperature at 
the transition to the thermodynamic limit and to find a quantitative estimate of this 
transition to the limit. Therefore, to complete the general picture, let us discuss the 
question of BEC as a phase transformation phenomenon. The order parameter here is 
the condensate density itbec = n /V : Ubec = for T > Tbec , and Ubec ~ 1 — T/Tbec 
for T < Tbec • The order of the phase transition is determined in a suitable classification 
according to whether a jump of the derivative of the heat capacity occurs upon transition 
of the temperature through T B ec • 

Expressions for the number of particles and energy are written without the surface 
terms as 



E = n e + V n k £ k = n £ + --3- P 5/2 ( ^— I • (2.23) 

l \ T \n + 1/ 



N = n ° 

k 



It follows from the definition (11 .11) of the temperature Tbec that 

V N 



*t C(3/2) 



(2.24) 



Introducing a normalized temperature t = T/T BEC and taking (I2.24p into account, we 
rewrite (01221) and fT2T23|) in the form 



C(3/2) J/z V^o + l 



n 



N = ri, + ——P V 2 — 7T » ( 2 ' 25 ) 



E n e 3 Nt 5 / 2 



n 



Tbec T BEC + 2^^\^+TJ " ^ 

It is seen from (12.251) that at temperatures in the vicinity of Tbec it ~ 1) an d under the 
condition 3> 1 , the number of particles in the BEC is also large: no ~ iV 2//3 , both 
above and below T B ec j m particular, for t = 1 



n, 



3/2 




^l^[l-0(N 3 / 2 )]. (2.27) 



C(3/2) 



For the specific heat c and its derivative at constant volume V , we obtain from (I2.25P 
and fl2T26|) 

1 (dE 1 dE\ dc_dc 1 dc 
°~ NT BEC \dt + udn~ ) ' dt ~ dt + udn~ ' ( ' ^ 



where 



For iV > n > 1 



.._ g _ 1 , 1 A/2[n /(«o + l)] ^ 
dn N-n n (n + 1) ^3/2(^0/(^0 + 1)] ' 



"4 + c(il)"« 3,2 + ( iV ~ 4/3 )- < 2 - 3 °) 
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This last expression shows that the scaling variables^ of the problem are 

C(3/2) nT 



T 



N{t 



V 



'71 N 

The asymptotic expansions for (12.251) . f)2.28p . and f)2.29p take the form 

2 1-y 



3 y 1 / 3 



+ 0(N 



-l/3\ 



dc 
~dt 



1 - 
3 



< i 27 

(yiV)l/3^ + l + y 

1 + Ay 



0(N 



-2/3\ 



1-7- 



where 



15C(5/2) 



a 



71 



+ 0(N^ 3 ), 

3 C 3 (3/2) 

7 _ 



(2.31; 



(2.32) 



4 C(3/2) ' ' C(3/2) ' ' IOtt C(5/2) 

It is seen from f |2.32j) that the small parameter in these expressions is iV~ 1//3 . Conse- 
quently, the average of physical quantities tends toward the thermodynamic limit rather 
slowly — according to an OfV -1 / 3 ) law rather than 0(V _1 ) as is usually assumed. The 
scaling regime sets in even more slowly, as is illustrated in Fig. 2. 




Figure 2: Derivative of the heat capacity as a function of the (scaled) temperature for different 
numbers of particles in the system: N = 10 3 (1), N = 10 6 (2), N — ► oo (3). 



At arbitrarily low but fixed deviations of the temperature from Tbec ( t — 1 + St) the 
variable y f)2.3ip has substantially different dependence on iV for different signs of the 
deviation St : 

" = ^ <1 " * >0 ' 

y = -N' ,3 6t > 1, St< 0. 



3 By scaling behavior we mean that some function (e.g., of two variables) f(x, s) , starting at certain 
scales s > sq degenerates into a function of only one, "scaling" variable y : f(x±, Si) = f(x2, S2) = f(y) 1 
where y = y(x, s) . In the problem under consideration the scale is the size of a system s ~ L ~ N 1 ^ 3 . 
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Accordingly, from f 1 2 . 3 2 j) we obtain the following expression for the jump of the derivative 
of the specific heat Ad = c'(l - St) - d(l + St) at N -> oo , St -> : 

Ac' = ^c max7 = T |-C 2 (3/2), 

2 l07T 

which agrees with the known expression given in [3j. Returning to the ordinary, unnor- 
malized temperature, we can see that the function d for iV — > oo goes over to a 9 
function with a "tooth" formed by the minimum on the curve (Fig. 2). 

Based on the given calculations, we stress that for BEC as a phase transformation 
phenomenon it is also necessary to have not too low temperatures at a fixed density of 
the Bose gas, so as to have a large number of (quasi)particles in the system@ Therefore, 
for experimental study of BEC it is preferable to increase the density of excitations at a 
specified value of T, as is done in Ref. [TS] . 

3 Bose condensation in thin films with anisotropic 
spin-wave spectrum 

A study of magnons in ferromagnetic microfilms with dimensions [H] 

L x ~ 10/xm, L y = L z ~ 1 (3.1) 

in a magnetic field directed along the surface of the sample, H||z , was reported in Ref. 
18 (see Fig. 1). The spectrum of long-wavelength spin excitations in real ferromagnets is 
rather complex. The corresponding dispersion relation e(p) differs markedly from (12.51) . 
which describes free massive particles. The additional contributions to e(p) are gener- 
ated by an external magnetic field and the magnetic dipole component of the interaction 
between spins. In the next Section the influence of these contributions on the BEC pro- 
cess are analyzed in more detail. In the simplest approximation one may keep only the 
fact that a magnetic-field gap appears in the dispersion relation, so that 

= (3.2) 

where Eh = ^hv-n is the gap energy, its frequency practically coinciding with the fer- 
romagnetic resonance (FMR) frequency. The presence of a gap in the dispersion relation 
does not affect the thermodynamics of the system, since this is only a shift of the en- 
ergy (frequency) scale. As is easily seen from the definition (12.91) of n , the substitution 
e — ► e + £h and \i — > fi + en leaves no , the number of particles on the lowest level, which 
essentially corresponds to the FMR, practically unchanged. And because it is no , not 
the chemical potential, that we use as the independent variable, all of the thermodynamic 
expressions of the previous Section remain in force. The quantity £h is conveniently used 
as a scale for parameters with dimensions of energy, if only because the diagnostics of the 
spectral density of of magnons in YIG in [18] was carried out at frequencies i/jj — 2 GHz , 

4 In this connection we note that an atomic Bose condensate with a total number of particles ~ 10 3 
(Refs. J3j and [5]), while qualitatively corresponding to the BEC phenomenon, nevertheless has little in 
common with the original and widely accepted conception of this phenomenon. 
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in the vicinity of the FMR. Then, e.g., the temperature is given as T = t/eu , which im- 
plies that, on the scale of quantities corresponding to those typical of [TS], a temperature 
of 1 corresponds to a value t ~ 10 . In the dimensionless quantities dispersion relation 
(13.21) takes the form 

e(p) = e k = e H (l + tOk) , uj k = '^2uj j k^ , (3.3) 

3 

where Uj = Ej/en, the Ej being defined in (12.51) . In the case of film dimensions (13. ip 
and magnon masses m m ~ 5m e , corresponding to YIG |10| . the numerical values of the 
parameters Uj are 

oj x ~ 10~ 4 , Lu y = u z ~ 10~ 10 . (3.4) 

In the experiments of [IB] , direct measurements were made of the magnon spectral 
density 

n(u) = ^11^5(00 — c^k) = ~ ; rrrr— — n (3-5) 

y 1 ~ e( w - w o)/*(l + no ) - 1 

with high resolution in the long-wavelength part of the spectrum, which is the most 
important and informative in the context of BEC. The magnon energy was determined 
from the frequency v = i>n(l + u) , and the number of magnons from the intensity of 
the emission in processes of inelastic Raman scattering of a light wave on the magnon 
distribution established in the film. 

On the right-hand side of (I3.5P we have introduced the spectral density of states 

g(u) = ^5(uj - uj k ). (3.6) 

k 

Because of the large difference in the longitudinal and transverse dimensions of the film, 
the ratio of the corresponding parameters ojj is also large: oJ x /uj y = uj x /oj z ~ 10 6 . 
This leads to the circumstance that the spectrum of states is split into layers around 
the harmonics corresponding to the first component of the quasimomentum. This pe- 
culiar structure of the spectrum allows one to decompose g(uj) and, hence, n(u) into 
three characteristic terms with different singular behavior in the low-frequency (near z/h ) 
region: 

g(w) = gc(w) + gi(u) + gaofa) , (3.7) 

where 

g c {uj) = 5{u -w ) (3.8) 
is the term corresponding to the contribution in (13.61) from the lowest magnon state; 



-u x - u y kl) , k 2 x = k 2 y + k 2 z , (3.9) 



is the contribution of the first excited layer (the fundamental harmonic of the first com- 
ponent of the quasimomentum) and, finally, 

oo 

g<x{u) = } j 5{u - u x kl - uj y k 2 ± ) (3.10) 

k x =2 k 
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is the contribution of all the other states. The prime on the summation sign in (13. 9p 
means that the term with kj_ = 2 is excluded, as it is taken into account explicitly in 
(13.81) . Changing from sums to integrals in (13.91) and (I3.10p . we obtain 



ex. 



TV f 7T 

gi{uj) ~ — / dkk 5{uj — u x — uj y k 2 ) = - — 6{oj — oj x — xuj y ) , (3-H) 

2 J ^UJy 



J y 

oo 

oo 



9oo{u) ~ 2^ dkkd(u -u x kl~u y k 2 ) = — ^6(0; -uj x kl) 





7T / I UJ 



AUJy V V UJ,, 





l)e{u-4w x ), (3.12) 



where < x < 3 , and the curly brackets on the right-hand side of (13.121) denote the 
fractional part of the quantity enclosed in them (i.e., {s} is the fractional part of the real 
variable s). Since < {s} < 1 , one can set the mean value of {s} ~ 1/2 , which gives 
for (I3TT2D 

g 00 {u) = - — ^— ( y/ui- -^[lo~ x J 9(u - Au x ) . (3.13) 

4ijJyy/uJx \ 2 J 

The partial contribution of g\{oS) (13. lip is proportional to the surface area of the sample, 
while in (00) (I3.13P the first term is proportional to the volume and the second to the 
area. Collecting the corresponding terms, we write g{uS) in the following form: 

g(u) = 9c{u) +gs v (u) + g v (u) , (3.14) 



where 



g c {u) = S(u-Uo), (3.15) 

7T 3 

9s v {u) = - — [0{uj-ujQ-HUJ y )--B{uj-Au x )\, (3.16) 

4tUy Z 



gvip) = -p^= 9{oj - Au x ) . (3.17) 

It is easily checked that for large u the phase volume element 

= [flv(w) +g Sv (uj)}duj 

agrees with the asymptotic estimate (12.141 However, in the low-energy (threshold ) region 
of interest to us, the question of the lower limits of integration over the phase volume 
becomes fundamental, and the answer to this question is contained in expressions ( 13.151) - 

(EUl). 

In sum, by substituting the formulas obtained for g{uS) into (13. 5p . we get 

n(u) = n c (uj) + n Sv {co) + n v (u) , (3.18) 
where, in accordance with decomposition (I3.14p . the functions 

n c {uj) = n 5(uj-LJo), (3.19) 

V{ ' 4uy (1 + n^ 1 ) exp(^) - 1 1 ' 

i \ * y/u 6(uj - 4uj x ) 

n v {uj) = — -, tt ; , (3.21) 

VK ' 4u yy /U£ (l + n - 1 ) eX p(^)-l V ; 
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specify the partial frequency distributions corresponding to the condensate, surface, and 
volume contributions. In a restricted range of frequencies u not exceeding a certain value 
uJ m ax , expressions (13.201) and (I3.2ip simplify if n 3> 1 and t ^> u) max : 

nt 9(u - u Q - >eu y ) - (3/2) 9(u - 4u x ) . . 

n S v {U) - — , [6.22) 

Auj y uo — u; + t/no 

riyfo;) = — —. — . (3.23) 

Aujyy/uJx oo — Uq + t/n Q 



It is seen from these expressions that because of the small parameter y/oJ^ in the de- 
nominator of (13.231) . the contribution (13.221) from the surface term in the general case is 
negligible: 

nv(uj) 3> ns v (oj) . But at the thresholds (u = ojo + xu y and uj = Auj x for (I3.22j) and 
(I3.23p . respectively) the situation is the opposite: ris v 3> riy . 

In real experiments the results of measurements depend on the resolving power of the 
device, and usually one observes not n(uj) but the spectral density averaged over some 
frequency interval 8uj ~ Q^ 1 : 

n obs {uj) = I du'n(u')F(u-J), lduoF{uj) = l ) (3.24) 



where F(uj) is the amplitudefrequency characteristic (AFC) of the receiver. Because 
the value of nc(oj) in (13 . 19[) is proportional to a 5 function, the observable condensate 
contribution 

n°g s (u) = n F(uj - u ) (3.25) 

actually reproduces the AFC. Expressions (I3.24p and (I3.25P show that the formation of 
the resonance peak as the BEC signal involves the participation not of just one state 
but also of other states whose energies are close to the gap £h • Therefore, at some 
distance from the "condensate" frequency ujq and not only at the thresholds, the partial 
contributions n°^(uj) and ny s (u) can have values comparable to n°^ s (uj) , despite the 
quantitative difference in the coefficients in expressions (13.221) and (I3.23p . Furthermore, 
owing to the finite width of the AFC the observable magnon spectral density n ohs {uj) is 
nonzero even in the region below the threshold ( to < , v < z/h ) • As the AFC we 
consider the often-encountered functional 

where Q is the quality factor. The integrals (I3.24p of the functions (I3.22p . (13.231) . and the 
Lorentzian (13.261) are taken explicitly. Taking into account the smallness of the parameters 



3 In radiophysics, for example, it corresponds to a circuit consisting of a single oscillatory loop. 
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(3.27) 



(3.28) 



(3.29) 



(3.30) 



where A denotes the uj -independent quantity (which depends on T , no , and the linear 
dimensions of the system) 



It can be seen from expressions H3. 28H — ( T3. 3 If) that the differential components of the 
spectral density behave differently upon variation of the temperature, film thickness, 
and the pumping done to produce magnons in sufficient number for BEC in the system. 
This satisfies the basic prerequisites for reliable separation of the condensate (coherent) 
component and the surface and volume (incoherent) components of the observed total 
(and, in essence, unified) spectral curve. 

It is extremely significant that in its main details the shape of the spectrum does not 
depend separately on the temperature and pumping but on their ratio t/no . Expressions 
fl3.28p - fl3.3ip show that the shape of the n obs (uj) curve is invariant with respect to simul- 
taneous variation of temperature and the number of particles in the condensate, provided 
that their ratio n = t/no is constant: at a fixed value of n variation of the temperature 
or magnon pumping affects only a coefficient that is common to all contributions. 

The critical magnon density Ubec (the derivative of the specific heat as a function of 
n ) jumps at the transition through Ubec ) has the form 



At room temperature T = 300 K and m m = 5 m e the thermal length (see (12.211) is small: 
At — 1-92 x 10~ 7 . Therefore the number of particles in the condensate is already quite 
large: ( N B ec/V — 6.28 x 10 14 , and the parameter r/ B EC = t/N B EC — 2.9 x 10~ 8 is 
sufficiently small as to ensure the validity of the approximations made in the derivation 
of fl3.28p - fl3.3ip . Despite the large number of particles in the condensate, its contribution 
to n obs (uj) is practically unnoticeable (roughly three or four orders of magnitude smaller 
than the thermal excitations). Nevertheless, please note that the n obs (u) curve already 
shows a completely formed peak as a consequence of the singular (at rj , uj —>■ ) behavior 
of the surface f)3.22p and volume (I3.23P contributions to n obs {uj) . For direct observation 
of n°(j s {uj) in the spectral density it is necessary that no 3> ./Vbec • 

It can be determined from (13. 280 - 03.311) that at such values of the parameter rj = r\ cros 
the number of particles in the condensate becomes macroscopic — e.g., equal to the volume 
contribution of the thermalized excitations. This gives a rough but simple estimate: 



A = 3 ln(t/no + 3u x ) — 2 \n(t/n + xuj y ) + In Q . 



(3.31) 



riBEc — 



C(3/2) 



(3.32) 




Alu, 



(3.33) 
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Figure 3: Spectral density of magnons and its partial contributions at different values of the 
parameter ry ; n ma x is the normalized coefficient corresponding to the maximum of the total 
spectral density at the given rj ; L x = 10 /xm , ryBEC — 4.8 x 10~ 9 , rj cros ~ 1.92 x 1CT 12 ; the 
symbols C, Sy , and V denote the corresponding partial contributions. 



For a film thickness L x = 10/jm this gives r] cros ~ 2.56 x 10 -12 . When the parameter 77 
reaches the value i] cros (with increasing magnon pumping) a kind of crossover occurs, i.e., 
the peak on the n ohs (uj) curve rises up sharply. A more accurate estimate for t] cros can 
be obtained from the quantity n^ s (0) = n° x s (u max ) , where n°^{uj) = ny s {uj) + n°^{ui) 
and LU max is the frequency at which the function n°^{uj) reaches a maximum. Hence we 
obtain r] cros ~ 1.92 x 10~ 12 for L x = 10/xm , and r] cros ~ 5.26 x 1CT 12 for L x = lfxm . 

^From these estimates one is convinced that at room temperature, t = 3 • 10 3 , this 
crossover occurs at no — 1.56 x 10 15 or a density no/V ~ 1.56 x 10 18 cm -3 , a value which 
has been attained in YIG films |18j . 

The above-discussed details of the behavior of the magnon spectral density as a func- 
tion of n Q (which corresponds to the dependence on the pumping level) are illustrated 
in Figs. 3-5. Figure 3 shows the spectrum corresponding to a microfilm (investigated in 
|18j ) of thickness ~ 10/zm, and Fig. 4 to a film an order of magnitude thinner, w lfirn. 
It is seen that in the thicker film at all pumping levels the surface contribution is less 
than the volume contribution, and, essentially, one can speak of only the competition of 
the condensate and volume contributions. The latter of these, being asymmetric with 
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respect to the maximum and rather strong in the region of the short- wavelength wing of 
the spectrum, weakens, yielding its place to the condensate contribution that ultimately 
determines the shape of the band at the highest pumping level (Fig. 3d). 

More interesting, perhaps, is the picture of the formation of the spectrum in a very 
thin film (Fig. 4), when all three contributions are present and competing. Moreover, 
the surface contribution, as we see, even at relatively weak pumping (Fig. 4a,b) not 
only remains larger but is essentially dominant, and together with the rising condensate 
contribution it largely determines the observed shape of the curve. Unlike the case shown 
in Fig. 3, here the volume contribution, which is also asymmetric, is relatively small 
(especially at large n ), and the condensate contribution mainly "combats" the surface 
contribution, prevailing only at the highest pumping levels (Fig. 4d). Figures 3 and 4 
clearly attest to the possibilities opened up in the study of BEC not only in thin magnetic 
films but also in other systems of finite size and with a small number of particles. 

Finally, Fig. 5 shows the evolution of the total distribution as the temperature is 
varied, leading, under otherwise equal conditions, to substantial growth of the number 
of magnons in the condensate with decreasing temperature. It should be kept in mind 
that the crossover point r] cros here cannot be interpreted as a phase transition signal. 
In the dimensionless variables used, the temperature of the latter, (11.11) . is given by the 
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expression 



fC_B-t BEC 4 1/3 
t-BEC — — — \WxWyUJ z ) 



AT 



2/3 

(3.34) 



.C(3/2) 

As was shown above (see Eq. (I3.32p ). the number of particles of the condensate at the 
phase transition point is A^bec ~ A^ 2//3 , and so 



^■BEC = £ (0J x UJyUJ z ) 1/3 

VBEC n BEC ~n C 2/3 (3/2) 
Vcros = C 2/3 (3/2 
?7bec tt 



2Q(u a p y w z ) 1/6 < 1. (3.35) 



This, in turn, indicates that the appearance of a noticeable peak in the spectral density 
essentially occurs already in the BEC phase. For the present case of a film with L x ~ 
10/im , L y = L z ~ 1 cm and for Q ~ 20 the ratio t] cros /rfQ^Q = N BEC /n cros ~ 4- 10~ 4 . In 
other words, at a fixed temperature a significantly lower density is necessary for a phase 
transition to the state with the condensate than for observation of the crossover, i.e., the 
transition to the state with a predominant number of particles in the condensate (see 
Figs. 3 and 4). 

It should be emphasized that the contribution n°gHoj) f 1 3 . 3 j) to the spectral density 
is a mesoscopic or finite-size effect that is isolated in analytical form. One notices that 
the frequency dependence of this contribution is practically no different from 
Meanwhile, the physical nature of these contributions is different: n°Q S {uj) is just the 
contribution of a large number of magnons with identical quantum numbers. This set of 
quasiparticles is a quantum object consisting of a macroscopic number of particles, or a 
coherent Bose condensate. On the contrary, ng*(oS) is formed of a set of magnons with 
different quantum numbers and corresponds to an obviously incoherent state. 

It is scarcely possible to distinguish these contributions by the attribute of coherence 
with the aid of the corresponding interference measurements. Nevertheless, because the 
contribution n°g s (u) is proportional to the temperature and the surface area of the sample 
[unlike n^ s (uj)}, these components of the spectral density can be distinguished with the 
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aid of a series of measurements at different temperatures on samples of different size and 
shape. 

For a quantitative description of the data it is necessary, among other things, to have 
information on the effective AFC of the measuring setup. The Lorentzian (13.261) that we 
used as a model can scarcely correspond to reality. Nevertheless, even in this simple model 
one sees a substantial dependence of the results of observation on the parameters of the 
AFC, in the present case — on the quality factor Q . Considering the partial contributions 
as functions of the reduced frequency Q uj , it is easy to see that the contribution of 
the condensate goes as ~ Q , while the volume contribution (13.291) goes as ~ Q 1 ^ 2 . 
Therefore, with increasing resolving power the contribution of the condensate becomes 
more and more noticeable against the background of the volume component. This can 
be used for identification of the volume component. However, the surface contribution 
is also practically proportional to the factor, and therefore variation of the Q factor is 
insufficient for separating it from the condensate component, and special data processing 
is necessary (such as that mentioned above). 



4 Inclusion of anisotropy of the magnon dispersion 
relation 

As we have said, in real ferromagnetic films, because of the contribution of the magnetic 
dipole interaction, the spectrum of long-wavelength magnons is anisotropic. Here the 
anisotropy is determined by the external magnetic field, which sets the magnetization 
direction [TD1 US] • Taking that interaction into account leads to a distortion of the isotropic 
quadratic dispersion relation (13.31) such that its dependence on the third component of the 
quasimomentum exhibits a characteristic dip (see Fig. 6) at a value p z = po — 3- 10~ 4 cm™ 1 
(recall that H||z). The corresponding e(p) dependence is well approximated by the 
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Figure 6: Form of the dispersion relation of long-wavelength magnons with wave vectors k||H 
in a ferromagnetic microfilm. 



expression 
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which for p a — > oo and e(0) — > £h goes over to (13.21) . If one knows the energy at zero, 
e(0) , and the parameters of the extremum, p and e = e(po) > then the momentum p a 
on the right-hand side of (14.11) can be written in the form 

P ~ = ^ - I- (4-2) 
Pa e - p z /2m 

In accordance with the phenomenological curve in Fig. 6 we have: e(0)/27rh ~ 4 GHz, 
Eo/27rh ~ 2.3 GHz. However, for the sake of unification in the case considered above we 
assume that the energy at the minimum, Eq , corresponds to the gap energy eh of the 
previous Section (e /2nh ~ 2 GHz). Then Eq. (14.11) is written in dimensionless variables 

as 

e(p) = £ k = e (l + lou), u k = u x kl + ujykl + uj z f(k z ), 

/(*.) = (k z - k f (l + -^-) = k\ + 2k f° + ka)2 - k (3k + 2k.), (4.3) 

\ K z + K a J K z + K a 

where 

k ~ 10 4 , k a ~ 270. 

We note that the minimum of energy in (14. 3 j) corresponds to the state with quantum 
numbers kk = (1, 1, /co) . 

Calculation of the spectral density of the number of states 

k 

with the anisotropic dispersion relation (14 .3p is complicated considerably in comparison 
with the simple calculations in Sec. III. Nevertheless, the final expressions are quanti- 
tatively very close to those obtained for the isotropic case. Repeating the calculation 
scheme of the previous Section, we separate out from g{uj) the lowest state gc(u) , the 
contribution of the first ( (k x = 1) ) harmonic g\{ui) , and the remaining part g^oS) : 

gfo) = gc{u) + g\(w) + fi'ocM, 
gc{u) = 5{u - uu x - u y ), 

gx{u) = ^ S[uj -u x - Uyk 2 y - uj z f(k z )} ~ 



J{u) 



6{u -u x - uj y ) . (4.4) 



The notation J(uj) on the right-hand side of (14.41) stands for the integral (which is eval- 
uated in Appendix B) 

m = [ , dz ( - s , (4.5) 

J sju) - uj x - w{z) 
where Z\ and z 2 are roots of the equation w(z) = u — u x : 

w{z) = u z f(z/^) = {z- z ) 2 (l + -^-) , (4.6) 
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where z a = k a ^fuT z ~ 2.7 • 10~ 3 , z = k ^/uT z ~ 1CT 1 . 

In an analogous way, by calculating the contribution g 0O (f jJ ') to accuracy O(^fuJ^) , we 
find 



#00 M 



Yl Yl ~ Uxk * ~ u y k l ~ Vtfi 1 **)] 

k x =2 k ± 

6{ljJ — Au) x — UJy) 



^(z 2 - zt) - -y/uJ^J(uj) 



(4.7) 



Collecting terms proportional to the surface and volume in (14.41) and (14.41) . we obtain (cf. 
< KW\> and (IXT7D ): 



- [9{w -u x -u} y )-- 9(u - 4u x - ujy)} , 

Z^UyUz 2 



7T (z 2 ~ Zi) 

V 



LO x LO y LO z 



(UJ - AUJ X - UJy) . 



(4.8) 
(4.9) 



It is easily checked that the integral J(oo) , generally speaking, depends weakly on u ; 
for example, for u ^> uj x the function J(oo) ~ vr/2, which coincides with our previous 
investigation of the isotropic case. A quantitative difference, though slight, appears only 
in the region of very low frequencies (near the threshold), as is illustrated in Fig. 7. 



1.15 




1.05 



Figure 7: Ratio of spectral densities of states: g a niyj) corresponds to the anisotropic dispersion 
relation (|4.ip . and gi s {uj) to the isotropic one ()3.2[) : as above, the symbol V denotes the volume 
contributions, Sy the surface contributions. 



Therefore, both the surface and volume contributions to the spectral density of thermal 
spin excitations 04. 7p are quantitatively similar from the case of an isotropic dispersion 
relation ( 13.2ft . On the whole, the observed spectral distribution in this case is qualitatively 
the same as the curves shown in Figs. 3 and 4. 

The calculations done lead to a completely optimistic conclusion for experiment: it 
would seem that the rather substantial difference in the dispersion relations (13.21) and ( 14.11) 
has almost no effect on the results of measurement of the spectral density of magnons, 
as is illustrated in Fig. 8. Thus the model of an ideal degenerate Bose gas with isotropic 
quadratic spectrum is completely adequate for interpretation of experiments on the BEC 
of magnons in real ferromagnetic systems. 
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Figure 8: The volume (a) and surface (b) contributions to the spectral density of magnons for 
a film 10 [im thick; the solid curve corresponds to the anisotropic dispersion relation (|4.ip . and 
the dashed curve to the isotropic one f|3.2|) . 



5 Conclusion 

The calculations done here convincingly demonstrate that in a thin ferromagnetic film with 
dimensions of 1 cm x 1 cm x 10 fim the BEC of spin excitations can be achieved already at 
a level of pumping that ensures the appearance and existence of ~ 10 15 magnons in the 
film. The fundamental thing is that there are practically no restrictions on temperature, 
and all the signs of BEC appear at rather high temperatures, including room temperature. 
Of course, lowering the temperature to the liquid hydrogen or helium region should lower 
the level of the "critical" pumping in view of the fact that, as we have shown, there is 
scaling or, in other words, the magnon spectral density depends not on the total number 
of particles and temperature but only on their ratio N/T . The shape of the observable 
spectrum of long- wavelength magnons remains unchanged when this ratio is held constant. 
We recall that the number of these quasiparticles is completely specified by the external 
pumping. 

On the other hand, the shape of the measured spectrum is largely connected to the 
AFC of the filter and can also be varied by improving its Q . Its increase (or the segrega- 
tion of a narrower frequency region around the lowest condensed magnon states) also leads 
to narrowing of the bandwidth, and consequently, to a more definite investigation specif- 
ically of the Bose condensate on whose consituent quasiparticles the Brillouin scattering 
of the light wave occurs. 

In interpreting the experimental data it should be recalled that one cannot claim direct 
observation of the condensate solely on the basis that a peak appears in the spectral 
density. All of the contributions have a similar structure, and it is necessary to separate 
each of them reliably from the experimental data. It follows from the theory that the main 
distinguishing features of the condensate contribution is its independence of temperature 
and the size of the sample. Therefore a precise and convincing separation of n £ S (uj) 
requires additional measurements at different temperatures. 
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In studying the BEC phenomenon in magnetic films, we have treated them as systems 
of finite size. In such systems, as we have said, it is known that BEC occurs not as 
a phase transition but as crossover. However, the condensation, which occurs near a 
certain temperature T cros , occurs so rapidly that that temperature, though not a critical 
temperature, can be treated as such. Apparently, the transition could be made more 
smooth (or smeared in temperature) if an even thinner film were used while the rather 
high measurement temperature was maintained. Here the absence of a phase transition 
point in no way precludes the phenomenon of BEC, or the macroscopic accumulation of 
magnons in their lowest energy state. Moreover, thin films are interesting in still another 
respect: in them one could trace the role and contribution of the correction terms (see Eq. 
(2.14)) arising when one goes from summation to integration over phase volume elements, 
something that, as far as we know, has hitherto escaped attention. 

The above calculations of the spectral density of magnons in a ferromagnetic film 
presuppose the presence of an equilibrium magnon gas, which, strictly speaking, does 
not correspond to reality. Intense electromagnetic pumping creates relatively short-lived 
spin excitations, as a result of four- magnon interaction, which, as was shown in [23] . 
rapidly relax to a quasiequilibrium distribution with a temperature equal or nearly equal 
to the temperature of the crystal (owing to the spinlattice coupling). Consideration of 
magnonmagnon and magnonphonon relaxation processes simultaneously with the process 
of magnon Bose condensation undoubtedly requires a special analysis and will be done 
separately. 

Finally, we note that it has been proposed in an experimental paper [18] that the 
accumulation of magnons at two (symmetric) points of the spin-wave spectrum fulfills 
the prerequisites for a nonuniform Bose condensate. However, from the standpoint of 
the thermodynamics of the process, no finite number of degenerate points of k space 
will affect the observed scattering pattern and the spectral density corresponding to it, 
which will be completely described in the framework of the approach stated above. This 
conclusion remains in force, despite the fact that the condensate that arises in such a case 
can be defined as incoherent. A distinction would arise only in the case when the lowest 
state of the magnon spectrum for some reason was infinitely degenerate (e.g., this could 
correspond to a model spectrum in the form of a trough). The condensate corresponding 
to that case is also incoherent, and the transition to it would have some distinctions. 
However, the study of the features of such a model is beyond the scope of this paper. 
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To demonstrate the problem with not separating the singular term, we consider a 
series whose sum is expressed in terms of elementary functions: 



APPENDIX A 




(A.l) 
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A "naive" transition from the sum to an integral in Eq. (A.l) gives the following approx- 
imation for S : 

oo 

S « S = [ = . (A.2) 

J P + x 2^ y 1 

o 

If the first term in Eq. (A.l), corresponding to I = , is separated off, and the remaining 
sum is approximated by an integral, we then get 



oo 



S~S l = 1 -+ / ^ = - + arctangyi) 
x J r + x x ^Jx 

1/2 

It is easy to check that approximation (A. 3) is practically no different from the exact 
expression (A.l) in the whole range of variation of the parameter x > . Approximation 
(A.2), on the contrary, is a very crude approximation compared to (A. 3), especially at 
x — > . For example, one can compare their asymptotic behavior at large and small x : 

1 7T 2 7T 4 X „, 0n 
^0: S = -+ - +0(x 2 , 
x o 90 

Si = - + 2-^ + 0(x 2 ), 
x 3 

Figure 9 shows graphs of S as a function of the parameter x for the exact and approxi- 
mate expressions. 




Figure 9: Behavior of the series S and its approximations So and S\ as functions of the 
parameter x ; the inset shows a plot of the ratio S/ Si . 
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APPENDIX B 



Consider the integral (see (14.51) ) 

m=\- r JtL= (b.i) 

J \/UJ — w(z) 
z\ 

where 

ro( , H(z _ 2o)2 ( 1 + _^_). 

The limits of integration in (B.I) are the positive roots (see Fig. 10) of the equation 

w(z lj2 ) =w- 




Figure 10: Graphical solution of the equation w{z\^) = oj . 



For calculating the integral, the expression under the square root in (B.I) is conve- 
niently written in the form 

(z 2 - z){z - z 1 )(z + z a + z 1 + z 2 ) 

u> — w(z) = . 

z + z a 

Then with the standard change of the integration variable 

_ zi(z 2 + z a ) + z a (z 2 - zi) sin 2 <p 
z — 



[z 2 + z a ) - (z 2 -^i)sin 2 v9 



integral (B.I) is transformed to 



tt/2 

2(z a + zi) f dip 



where 



J(u) = . v a 17 / y . , (B.2) 

V (z a + z 2 )(z a + 1z\ + z 2 ) J (1 - asin 2 ^)^! - absm 2 (p 



z 2 — Z\ Z 2 + Z\ 

a = , o 



Z a + Z 2 Z a + 2Zl + Z 2 
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The integral on the right-hand side of (B.2) is the Legendre complete elliptic integral of 
the third kind, II (a, ab) . 

NOTE ADDED IN PROOF 

After this article was submitted to Low Temperature Physics, an experimental paper 
appeared [J. Kasprzak, M. Richard, S. Kundermann, A. Baas, J. M. J. Keeling, F. M. 
Marchetti, M. H. Szymanska, R. Andre, J. L. Staehli, V. Savona, P. B. Littlewood, B. 
Deveaud, and Le Si Dang, Nature 443, 409 (2006)] in which the BoseEinstein condensation 
of quasiparticles was also reported. In that paper a different type of relatively light 
quasiparticles of the Bose type was considered, specifically, exciton polaritons, a high 
density of which can be produced by laser pumping in optical microcavities in a CdTe 
crystal. In it, as in YIG, thermalization of the pumped quasiparticles occurs, and above 
a certain quasiparticle density there is macroscopic occupation of the quasiparticle (in 
this case, polariton) ground state, at a rather high critical temperature Tbec ~ 19 K. 
Evidence of Bose condensation is provided by weak coherent correlation effects observed in 
the radiation emitted from the crystal. There, however, it is not ruled out that additional 
information about the formation of a Bose condensate specifically could be obtained if 
the volume and surface contributions to the observed intensity and line shape could be 
separated (on the basis of their temperature or concentration behavior). Although in our 
paper these contributions are calculated for magnons, they should undoubtedly also be 
present for polaritons in the small volumes of the optical microcavities. 
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